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Abstract 

The Hamiltonian Hopf bifurcation has an integrable normal form that describes 
the passage of the eigenvalues of an equilibrium through the 1 : — 1 resonance. At 
the bifurcation the pure imaginary eigenvalues of the elliptic equilibrium turn into a 
complex quadruplet of eigenvalues and the equilibrium becomes a linearly unstable 
focus-focus point. We explicitly calculate the frequency map of the integrable normal 
form, in particular we obtain the rotation number as a function on the image of the 
energy-momentum map in the case where the fibres are compact. We prove that 
the isoenergetic non-degeneracy condition of the KAM theorem is violated on a curve 
passing through the focus- focus point in the image of the energy- momentum map. This 
is equivalent to the vanishing of twist in a Poincare map for each energy near that of 
the focus- focus point. In addition we show that in a family of periodic orbits (the non- 
linear normal modes) the twist also vanishes. These results imply the existence of all 
the unusual dynamical phenomena associated to non-twist maps near the Hamiltonian 
Hopf bifurcation. 

Keywords: Hamiltonian Hopf Bifurcation; KAM; isoenergetic non-degeneracy; Van- 
ishing Twist; Elliptic Integrals 

1 Introduction 

Understanding the dynamics of a Hamiltonian system near equilibrium points is of fun- 
damental importance. In the elliptic case the eigenvalues of the linearization are pure 
imaginary, Xj = ±iujj, j = 1, ... , n, where n is the number of degrees of freedom, which 
will be 2 in the following. By the Lyapunov Centre theorem the n normal modes of 
the linear approximation persist in the non-linear system when the eigenvalues are non- 
resonant. The resonant cases were much more recently treated by |2L)1 1151 El EI- In 
some sense the most exceptional resonance is the so called 1 : — 1 resonance, in which the 
quadratic part H2 of the Hamiltonian H has degenerate eigenvalues and is not definite, 
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2H 2 = uj(p 2 +x 2 ) — u>(p 2 +y 2 ). The unfolding of the normal form gives a family of Hamilto- 
nian systems with an equilibrium point that looses stability by passing through the 1 : — 1 
resonance. This is called the Hamiltonian Hopf bifurcation, and was studied in detail by 
|18l I19j . also see (21 EI- The corresponding family of Hamiltonians is Liouville integrable. 
In this paper we show that near this bifurcation the twist vanishes. This means that the 
rotation number (i.e. the ratio of frequencies of the Hamiltonian flow on invariant tori) of 
the one parameter family of invariant tori with fixed energy has a critical point. 

The dynamical consequences of vanishing twist are well known. They were first de- 
scribed by |12| . and later studied by HI ITTj. In ^U] we have shown that the vanishing 
of twist in one parameter families of maps occurs near the 1 : 3 resonance, also see |13| . 
In P| we have shown that also in 4 dimensional symplectic maps the vanishing of twist 
appears near resonance. More recently in [S] we have shown that the twist also vanishes 
near the saddle-centre bifurcation, in which one multiplier passes through zero. In this 
paper we show that the principle that the twist vanishes near resonance also applies in 
the Hamiltonian Hopf bifurcation. For flows the condition of vanishing twist is one of the 
conditions for the standard form of the KAM theorem to hold. In this setting it is usually 
called the isoenergetic non-degeneracy condition. There exist KAM theorems with weaker 
conditions El , so that vanishing twist does not necessarily imply that the torus will 
be destroyed. It does mean, however, that a resonant twistless torus will create all the 
unexpected dynamics described by the twistless standard map. 

In the following two section we present well known material about the Hamiltonian 
Hopf bifurcation, in order to introduce the Hamiltonian and its Energy-Momentum map 
and to fix our notation. Then our own contribution starts with the derivation of the 
actions and the rotation number. The rotation number and its derivative are analysed 
near critical values of the energy-momentum map, namely near the isolated focus-focus 
point in the compact case and on the family of relative equilibria. The details of the 
expansion of the elliptic integrals are given in the appendix. 

2 Hopf Normal Form 

Consider coordinates q = (qi,q 2 ) an d conjugate momenta p = (pi,p 2 ) so that the sym- 
plectic form on M 4 is Q = dpi A dqi + dp2 A dq 2 . The normal form for the Hamiltonian 
Hopf bifurcation is 

H( qi ,q 2 , Pl ,p 2 ) = pr x + r 2 + 5( 7 ri + r 3 ) + erf + 2Bt 1 t z + 2DT 2 + o 3 (r t , r 2 , r 3 ) , (1) 

where T% = p 2 q\ ~Piq2, ^2 = ^3 = 5(^1 + an d 5 is a bifurcation parameter, 

/3,7, B,C, D are real constants such that 0^=0 and D / 0. The expression O3 denotes 
terms of order no less than 3 with respect to Ti,i = 1, 2, 3. For simplicity we will use the 
notation u = + 5j. The system has an equilibrium point at the origin Pi = qi = with 
eigenvalues \J — 5 ± iuj. For ease of notation we write a = y— 5 when 5 < 0, so that the 
eigenvalues of the equilibrium are ±a ± iu. The dependence of u> on 5 is not essential for 
our purposes, because 7^ 0. By a symplectic scaling with multiplier the parameter uj can 
be scaled to 1, and D can be scaled to ±1 at the same time. We find it useful, however, 
to keep unsealed variables and parameters until the very last section. 

The Hamiltonian system Q is Liouville integrable. A second independent constant of 
motion is T\. It generates the S 1 symmetry 

^^xl^l 4 , <f>(#,q,p) = (S M ,Sop), S*=(_°^ ^J) (2) 
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Figure 1: A typical intersection M^j of the surfaces Pj and H- (h). 

The corresponding momentum map J : R 4 — > R is given by J(q,p) = P2Q1 — P1Q2, which is 
T\. Since J generates the periodic flow $ with period 2ir it is an action of the integrable 
system. We denote the (constant) value of J by j. To perform the reduction with respect 
to this symmetry we use invariant theory, see e.g. Singular reduction occurs in this 
example because the action $ is not free: the equilibrium (= the origin) is a fixed point 
of this action. The algebra of polynomials in R 4 that are invariant under $ is generated 
by Ti, T2, T3 and T4 = piqi + p2(?2- This means that any polynomial of qi, q2, Pi, P2 that 
is invariant under $ can be written as a polynomial of Tj, i = 1, ... ,4. The generators 
satisfy the relations 

G(T) = rf/2 + T 2 J2 - 2T 2 T 3 = 0, T 2 > 0, T 3 > 0. (3) 

The reduced phase space Pj = J~ l {j)/ S 1 is defined by (J3J) with Ti = j as a semialgebraic 
variety in R 3 with coordinates (^,^^4). If j ^ the reduced phase space Pj is one 
sheet of a two-sheeted hyperboloid given by (3), so it is a smooth manifold. But for j = 
it is half of an elliptic cone and hence is not smooth because of the singular point of the 
cone at the origin (r2,T3,r4) = (0,0,0). 
The reduced Hamiltonian is 

Hj (r 2 , r 3 , r 4 ) = w j + r 2 + <5r 3 + cf + 2Bjr 3 + 2DV 2 3 . (4) 

The surface HJ l (h) is a parabolic cylinder in (r2,T3,r4) that is independent of T4. The 
integral curves of the reduced system are given by the intersection of the surface Pj with 
the surface HJ (h), as illustrated in Fig. ^ We denote the intersection by Mhj. 

The Poisson bracket {•, •} associated with the standard symplectic structure fi on R 4 
defines a Poisson structure on the algebra of invariant polynomials with brackets 

{r 1 ,r 2 } = {r 1 ,r 3 } = {r 1 ,r 4 } = o, (5) 
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{r 2 ,r 3 } = r 4 , {r 3 ,r 4 } = -2r 3 , {r 4 ,r 2 } = -2r 2 . (6) 

The momentum map of the S 1 action used for the reduction induces the Casimir Ti in 
this Poisson bracket. In addition also the relation between the generators © given by 
G is a Casimir. Accordingly the nonzero brackets give a Poisson structure on R 3 with 
coordinates (r2,r3,r4), that has the reduced spaces Pj as its symplectic leaves. It can be 
written as 

{r 2 ,r 3 } = dG/dr 4 , {r 3 ,r 4 } = dG/dr 2 , {r 4 ,r 2 } = dG/dr 3 . (7) 

The reduced equations of motion are 
T 2 = {r 2 ,i^} = r 4T^r-, 

f 3 = {r 3 ,^} = -r 4 , (8) 

dHj 



r 4 = {r 4 ,# i } = -2r 2 + 2r 3 



The integral curves of this flow are given by Mhj, the intersection of Pj and HJ 1 ^). 
In general the intersection of these two manifolds is either empty or diffeomorphic to a 
circle. The preimage of any point in reduced phase space is the set of points in original 
phase space that are mapped to this point by the momentum map J. If one point in the 
preimage is known the others can be obtained by letting the flow of J (i.e. the map <&) act 
on this point to get the complete fibre. This gives a circle unless starting in the origin. 
Therefore the preimage of a circle M^j is a two dimensional torus T 2 in the phase space 
of the original system. 

Exceptions occur for equilibrium points of the reduced system. They occur either when 
the surface HJ (h) is tangent to Pj or when j = and Hq (h) contains the singular point 
at the origin (which implies h = 0). The preimage of the singular point is not a circle, 
because T 2 = T 3 = implies q\ = g 2 = p\ = p 2 = and this is a fixed point of the flow 

This is the equilibrium point in the full system that undergoes the Hopf bifurcation. 
All other equilibrium points of the recuded system are reconstructed to periodic orbits of 
the full system; they are relative equilibria of J. The S 1 action generated by J is not free. 
The origin is a fixed point, and this is the reason why singular reduction is needed in this 
example. 



3 Energy-Momentum Map 

Using the reduced system we can find the critical values of the energy momentum map 

F:R^R 2 , (p,q)» F(p,q) = (H(p,q), J(p,q)) . (9) 

The values of the energy-momentum map are denoted by (h,j). For every regular value 
of F the preimage in phase space is a two dimensional torus. The critical values are 
determined from equilibrium points of the reduced system because their preimages are 
not T 2 . Since we are interested in a neighbourhood of the origin in phase space for small 5 
we will only consider a small neighbourhood of the origin in the image of the momentum 
map. 
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Consider the reduced equilibrium points caused by the singularity in the reduced space 
first. This singularity occurs for j = 0. The singular point (T2,r3,r4) = (0,0,0) has 
energy h = 0. The equilibrium at the origin in phase space is therefore mapped to the 
origin in the image of the momentum map. 

When 5 > the intersection Mq^ restricted to a neighborhood of the origin in reduced 
phase space consists only of the origin. It reconstructs to an elliptic equilibrium. However, 
if 5 < then Mo,o is a non-smooth circle with a corner, if it is compact. The preimage of 
Mo,o is diffeomorphic to a pinched torus in this case. 

Consider next the equilibrium points caused by a tangency of Pj and Hj l {h). At 
these critical values of Hj the gradient of Hj and the gradient of G are parallel. Since 
dG/dT^ = T4 the tangency may occur only on the hyperplane 1?4 = 0. The intersections 
of Pj and Hj l {h) with this hyperplane are one branch of a hyperbola and a parabola, 
respectively. They are described by the equations 

j 2 = 4r 2 r 3 , (10) 

h = ojj + Cj + T 2 — h 2D T 3 H — - — . (11) 



8D V 40 

At the extremal values of h the two curves are tangent. Eliminating r 2 in (|10j) using (|11|) 
gives a polynomial of degree 3 in T 3 depending on j and h 

Q 3 (T 3 ) := -8DTl - 4(5 + 2Bj)T 2 3 + A(h - uj - Cj 2 )T 3 - j 2 = . (12) 

This polynomial Q 3 (T 3 ) gives the value of T 2 obtained from G = and expressed in terms 
of T 3 . The tangency between the hyperbola (|10|) and the parabola (|1 1 1) occurs when Q 3 
has a double root. We will first discuss all values of (h,j) for which a tangency occurs, 
irrespective of them satisfying the constraints T 2 > and T 3 > 0. In a second step the 
critical values of the energy momentum will be found by consideration of these constraints. 

To parametrize all tangencies we make the ansatz Q 3 (z) = — 8D(z — d/2) 2 (z + s 2 / D /2) 
with parameters d and s parametrising the double and single root of Q 3 , respectively. This 
leads to the parametrisation of the tangencies by s S R 

icritOO = «*(-), d (*) = 2 *B~8-Dy (13) 
hcrit(s) = wjcrit(s) + Cj crit (s) 2 - 2d{s) 2 D + 2sj crit (s) • 

The root T 3 = —s 2 /D/2 always has the opposite sign than D. The curve (h CT it(s), j C rit(s)) 
has singular points when s has one of the singular values satisfying 2Bs 3 — 3s 2 D + 5D = 0. 
The number of singular points changes when the discriminant 1085 D 2 (D 2 — 5 B 2 ) vanishes. 
For small \5\ the only change occurs at 5 = 0, see Fig. [2] for two slices of the "swallowtail". 

For 5 > the curve has two singular points near the origin for some s E (—y/S,y/6). 
The two singular points are located at 

h* = ±^5 3 / 2 + WU -J D 5 2 ± 0(5 5 / 2 ) 

V27D 18D 2 (i4) 

j* = ±^^5 3 / 2 + -^8 2 ± 0{5^ 2 ) 
J V^ID 9D 2 V ; 

for small 5 > 0. The curve of critical values has a self-intersection at s 2 = 5. The 
intersection point at the origin marks the elliptic equilibrium with eigenvalues i(ui ± yS). 
The slopes of the intersecting curves are given by the imaginary parts of the eigenvalues. 
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Figure 2: The discriminant surface of the polynomial Q3. The critical values of the 
energy-momentum map are contained in the discriminant surface. The small triangular 
part for 6 > is not in the image of the energy-momentum map for D < 0. The figures 
correspond to the cases of S = —0.15 and 5 = 0.15 respectively. (B = C = —D = oj = 1.0) 



6 



For 5 < the equilibrium point at the origin is unstable with eigenvalues y/—5 + iu. 
The curve (^ C rit(s),icrit(s)) does not intersect the origin, instead the origin is now an 
isolated critical point. The curve is above the origin for D > and below for D < 0, e.g. 
the point with s = is at (h,j) = (-5 2 /8/D, 0). 

Only tangencies that occur with the part of the hyperbola in the positive quadrant 
give critical values of the energy-momentum map since T2 and T3 are both non-negative. 
The double roots d(s) near s = occurs at 

d(a) = ~-0(s). (15) 
Near the intersection at the origin around s = ±5 the double root is 

d(s) = ^(±«- V5) + 0((±a- v^) 2 ) (16) 

If D > this implies that the smooth curve for 5 < is in the (boundary of the) image of 
the energy-momentum map, while the part of the curve with s 6 (—yS, yo) is not in the 
image for 5 > 0. Conversely for D < the smooth curve for 5 < is not in the image, 
while for S > only the triangular part with s & [~VS, y/S] is in the image. Therefore the 
union of the bifurcation diagram (i.e. the set of critical values of the energy-momentum 
map) for +D and —D gives the discrimiant of the polynomial Q3. 

The type of the preimage of the critical values is determined by the character of the 
intersection Mhj- The positive half of the hyperboloid G projects onto the area above the 
hyperbola given by j 2 < 4T2r3. If the parabola (|llj) touches the boundary of the area 
from the outside, the preimage in the full phase space is a circle, hence a stable periodic 
orbit. If the parabola touches from the inside, the preimage is a circle with a separatrix, 
hence an unstable periodic orbit. This can only occur when D < 0, because then the 
parabola is open upwards. It only occurs when 5 > for s between the two singular values 
enclosing zero. In all other cases the parabola touches from the outside. The complete 
bifurcation scenario in the two cases therefore is as follows, see Figure |2] for illustration. 

The case D > 0: For 5 < there is an isolated focus- focus point at the origin and a 
smooth curve of elliptic periodic orbits nearby. For 5 > there is an elliptic equilibirum 
point and there are two families of elliptic periodic orbits (non-linear normal modes) 
emanating from the equilibrium. 

The case D < 0: For 5 < there is nothing but an isolated focus-focus point. For 5 > 
there is an elliptic equilibirum point and there are two families of elliptic periodic orbits 
(non-linear normal modes) emanating from the equilibrium. Both families terminate in a 
cusp formed with the same family of hyperbolic periodic orbits. The set of critical values 
therefore forms a triangle with two cuspoidal corners and one regular corner at the origin. 



4 Actions 

From now on we shall assume the parameter D to be a positive number. In this case each 
constant energy level is compact and by the Liouville-Arnold theorem it is possible to 
define action-angle coordinates near regular points of F. To construct the second action 
we need to integrate a 1-form O over Mjh, where the differential d0 coincide with a 
symplectic structure induced by the quotient map n^j : M 4 — > F~ 1 (h,j)/S 1 from the 
original phase space, i.e. d0 = {t:^)Vi\t M h j ■ To find the form we chose T3 as one 
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variable and find its conjugate variable by solving the equation 



{/(r 2 ,r3,r 4 ),r 3 } = i. (17) 

A solution of (|17|) is the function 

/(r 2 ,r 3) r 4 ) = ^. (is) 

ZL 3 

Then the canonical one form is O = 7^-dr 3 and we obtain 



r(h-\ 1 / r4dr3 1 I V 41 ^-; 2 ^ hQ , 

M h,j Mh,j 

for the second action. Here r 4 is considered as a function of r 3 by first expressing r 4 in 
terms of r 2 and T 3 on the reduced phase space Pj and then by expressing r 2 in terms of 
T 3 using Hj = h. As a result the polynomial Q3 is found as already given by (|12j) . The 
action integral hence is defined on the elliptic curve 

£ = {(w,z) :w 2 = Q 3 (z)}. (20) 

Recall that w = T 4 on Pj. Now the action integral can be written as 

J 2 = ^-{-dz. (21) 
An J z 

It is an integral of the third kind with a pole at z = and residue ±17/471". 

The formula for the action can also be obtained in a classical way, using polar coordi- 
nates as in |18j . A slightly different coordinate transformation illucidates the connection 
between the two approaches. The new symplectic structure is f2 = dP g A dg + d J A dcf) 
and "symplectic polar coordinates" valid for q\ + q\ > are introduced by 

/ — / — sin 1/? 

qi = y , 2gcosip, pi = P g y/2g cos ip - J— 7=, (22) 

V z 9 



1 — / — cosy? 
Q2 = y / 2gsmip, p 2 = P g \J2g sin ip + J—t=- (23) 

V z 9 



The invariant polynomials are related to these coordinates by 

Ti = J, T 2 =gP 2 + — , F 3 = g, T i = 2gP g . (24) 



In these variables the Hamiltonian takes the form 



72 

H = gP 2 + — + u J + (6 + 2BJ)g + CJ 2 + 2Dg 2 , (25) 
y 4# 

and the equations of motion are 
<p = -^-+uj + 2Bg + 2CJ , 

9 = 2gPg , (26) 
J 2 



P 9 = P 2 g -i^ + {5 + 2BJ)+ADg 
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Solving the Hamiltonian l|25jl for Pg gives 

P 2 = QM (27) 



so that the action integral ()19l) is obtained from integrating the canonical form P g dg over 
a path with constant (f). 

5 Rotation Number 

We want to check the isoenergetic non-degeneracy condition of the KAM theorem. A 
torus is non-degenerate in this sense if the map from the actions restricted to a constant 
energy surface H = h to the frequency ratios oj\ : 0J2 is non-degenerate. This means that 
the frequency ratio (or rotation number) W = uj\/u)2 changes when the torus is changed 
at constant energy. On a local transversal Poincare section this condition is called twist 
condition. 

In our case this is equivalent to the non-vanishing of the partial derivative of the 
rotation number W with respect to the action J\. By definition the winding number is 
the ratio of frequences lo\,oj2, corresponding to the actions J\ and J2. If the Hamiltonian 
is expressed in terms of J = J\ and J2 then d\H( Ji, J2) = u)\ and d2H(J\, J2) = ll>2- 
Therefore we find by implicit differentiation of J2(h,j) = J2 that 

W = — = --rrr- (28) 

However, the simplest way to obtain W is to observe that it is the advance of the angle 
4> conjugate to J during the time of a full period of the motion of g = IV The period of 
the motion is obtained from the reduced equation of motion I^ = — T^. On Pj this gives 

f| = 4r 2 r 3 -j 2 (29) 

and eliminating T2 by using Hj = h gives 



dlV 2 



4r 3 r 2 (r 3 Aj)-j 2 . (30) 



dt 

By separation of variables we obtain the period of the reduced motion as 

To obtain the advance of <j> in time T we change the time t in (|26|) to "time" T 3 and find 
d0 j + 2r 3 (cu + 2Ci)+ J B(2r 3 ) 2 



dr 3 2r 3 r 4 



(32) 



Expressing in terms of T 3 on the reduced phase space Pj as before the period of the 
solution of this equation gives the rotation number 2nW . The rotation number W can 
therefore be written as a linear combination of integrals of the first, second and third kind, 

2irW(h,j) = (uj + 2Cj) I — + 2B([ — + J -([ — . (33) 

w w 2 zw 
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The first integral is of the first kind and proportional to the period T. The last integral is 
of the third kind and propotional to the action J2. When D > the polynomial Q3 has 
three real roots, which we denote by z neg , z m - m , z max , such that z neg < < z min < 2 max . 
The closed loop integrals encircle the finite range of positive Q3, and therefore can be 
rewritten by the rule 

(34) 

The elliptic integrals can be transformed to Legendre standard integrals K(k), E(k), and 
H(k) of the first, second, and third kind, respectively, with modulus k and characteristic 
(or parameter) n given by 

jj2 ^max Z m [ n ^ ^max -2-min (35) 

Zm&x -2-neg Z mSiX 

The result is 

w ,. ^ = (2Bz neg + uj + 2Cj)K{k) + 2ff(z max - z neg )E(k) + jU(n, k)/(2z max ) 

7T\/2_D^/Z max Z neg 

Explicit formulas for the vanishing of the twist dW/dj can be derived from this. 

The rotation number is a complicated function of the constants of motion (h,j). The 
level lines of this function are shown in Fig.|3]and|IJ Near the cases where the discriminant 
of the elliptic curve E defined by Q3 in (|20|) vanishes, simpler formulas can be derived. 
This occurs either at the boundary of the image of the energy-momentum map described 
by (|13|) . or at the isolated focus-focus point (h,j) = (0,0) inside the image. In the next 
section we will treat the latter case. 



6 Rotation Number near the Focus-Focus Point 

We introduce a small parameter by scaling h and j by epsilon, hence replace h — * he, 
j — > je. This means that we obtain an expansion that approaches the origin on a ray. 
Alternatively one can view e as a formal expansion parameter that keeps track of the fact 
that both, h and j are small and of the same order. The focus point only exists for 5 < 0, 
which we henceforth assume. At the origin Q3 reduces to 

Q 3 (z;h = 0,j = 0) = -Az 2 (2Dz + 5), (37) 

so that z m j n and z ncg collide at and z max = —S/2/D. The roots can be expanded in 
power series in e, and the result is 

2z neg = - l -±P £ + 0(e 2 ), 
a 

2z min = - l -^e + 0(e 2 ), (38) 
a 

q 2 Dl — iBa „ . 9 . 
2z max = — + 2 -f- e + O (e 2 ) . 

Here and in the following we use the abbreviations 

; h ~ u J 2 ;2 , -2 f D r~x ron\ 

I = , p A = V + j, f = a = V-6 . (39) 

a a- 3 
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-0.0008 -0.0096 -0.0004 -0.0002 0.0002 0.0004 0.0006 0.000Q 

Figure 3: The level lines of the rotation number on the (J, h) plane. The bifurcation 
parameter 5 is decreasing from top to bottom, 6 = —0.2,-0.15,-0.02 (u = B = C = 
—D = 1). The level lines are spirals. The boundary of the image of the energy momentum 
map given by a family of relative equilibria is also shown. 
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Figure 4: The level lines of the rotation number on the (j, h) plane. Here the imaginary 
part of the eigenvalue ui is changed, from top to bottom uj = 1/2, 1,2 (5 = —0.15,5 = 
C = -D = 1). 
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( 4 °) 



The expressions are only real in the case 5 < 0, otherwise the focus-focus point does not 
exist. Inserting this into ()35|) gives 

k 2 = 1 - 2fpe + 0(e 2 ^ 
n = 1 + f(i - p) £ + 0(e 

For small e both, fc and n are close to 1 and they satisfy the inequality 

k 2 < n < 1 . (41) 

In the limit k — ► 1 the elliptic integrals are singular, but there are expansions that include 
the logarithmically diverging terms. The details of this expansion can be found in the 
appendix. The result for the rotation number is 

2nW(h,j) = --hip-tan" 1 ~ + O(l) . (42) 
a I 

Keeping terms only up to order 1 is enough because when the twist condition dW/dj is 
calculated the present terms both give singular contributions, the constant term disappears 
and the first order term in e is very small compared to the singular terms. Note that W is 
not a single valued function. The fact that tan _1 (j//)/27r changes by one when the origin 
is encircled is an expression of the monodromy of this focus-focus point, see [7j. The 
expression given by the elliptic integrals, see the appendix, gives a continuous function 
which is, however, not differ entiable when j = 0. 

The level lines of W are spirals, which are easily parametrized in polar coordinates 
with the radius p as parameter. Instead of viewing these spirals as the level lines of W 
where W is a many-valued function, it might be easier to see them as the integral curves of 
a flow in the plane that has an equilibrium point at h = j = 0. The linear approximation 
to this flow is best written in complex notation z = I + ij so that z = \z where A = a + iu>. 
Therefore we may formulate the result like this: The level lines of W are the integral 
curves of a planar node with the same eigenvalue as the focus-focus point. 

Calculating the leading order condition for the vanishing of twist is easily done by 
differentiating the first two terms of the expansion with respect to j. The result is that 

= h{u 2 - a 2 ) = ju(uj 2 + a 2 ) . (43) 

This means that the twist vanishes on a line that has a tangent whose slope at the focus 
point can be read of from the above expression as 

h uj 2 + a 2 , AA . 
~ = 2 • ( 44 ) 

In particular this slope only depends on the eigenvalue of the focus- focus point. Ap- 
proaching the bifurcation point a — > 0, and hence the slope approaches uj, see Fig. Eland 
® 

7 Vanishing Twist of Periodic Orbits 

Now we consider the vanishing of the twist for the relative equilibria (/i C rit(s) ; Jcrit(s)) in 
the boundary of the image of the energy-momentum map. First we rescale the parameters, 
in order to keep the formulas manageable: 

Db Dc 5 . (lu 2 

B = — , C=— , p = — , J = ~rr, h = -FT> Z = -FT- ( 45 ) 
U) uj z uj z D D D 



13 



Note that when 5 < 0, which is the case we are interested here, we have [i = —a 2 /to 2 . In 
terms of the new parameters the polynomial Q(z) becomes 

Q(z) = ^Q(C), Q(C) = -(C 3 + (A* + 26j)C 2 + 2(j + cj 2 - A)C + j 2 ), (46) 
and the derivative dW/dj is 

kl »(/ a+gwc + / (i»c 2 +(i+2ej)c+ J ) 2 d C \ {47) 



dj W c<3 1/2 (0 / C<3 3/2 (0 

The parametrisation of the critical values has already been obtained in (|13j) . After the 
scaling it reads 

d ( s ) = oru , TT > icrit(s) = sd(s), /icrit = w j crit (s) +cj 2 rit (s) +2(i(s) 2 + 2si crit (s) . (48) 

On the curve (j C rit(s), ^crit(s)) the integrals on the right hand side of (|4T|) can be computed 
by the method of residues. Replacing j and h wherever they appear by their critical values 
parametrised by s gives a condition for the vanishing of the twist on the curve of critical 
values. It implies the following polynomial equation 

9 

R 9 (s) = Y,a k s k = 0, (49) 

fc=0 

where 

a = - 12/i - 4(4 + % 2 - (3b 2 - 16c)/i 3 , 

oi = - 24(5 + b)p - 4(6 + b 2 - 6c)fi 2 + 46c/i 3 , 

a 2 = - 84 - 12(1 + 146 + 6 2 )/i + 3(156 2 - 8c + 86c)/x 2 - 12c 2 fj 2 , 

a 3 =4(30 - 666) + 16(6 + 6 2 + 9c) fx + (246 3 - 526c)/i 2 , 

a 4 =12(5 + 256 - 236 2 ) + (646 3 - 296 2 + 192c + 3366c)/; - 4c(166 2 + 15c)/i 2 
05 = - 12(86 3 - 316 2 - 76 + 14c) - 126(126 2 - 29c - 166c)/i - 966cV, 
a 6 =32 06 3 + 836 2 - 31 2c - 3606c - (966 4 - 646 2 c - 156c 2 
a 7 =46(326 3 + 386 2 - 171c - 486c) - 646c(26 2 - 3c) fi, 
a 8 =12(86 4 - 326 2 c-7c 2 ), 
ag = — 966c 2 . 

The first two coefficients ao and a\ are of the first order in fj,, while a 2 is non-zero for small 
fi. Therefore the roots of Rg can be expanded in powers of \fi\ 1 ^ 2 = a/ui, and for small \fj,\ 
there are only two roots in a l/ij 1 ^ 2 neighborhood of the origin. The result is 



(50) 



S± = ±-y/^ + 1(10 - 6> + 0{\^' 2 ). (51) 
The critical values of the twistless periodic orbit are then given by 



fc± = ±~]/-f + |(43/2 - 6% 2 + 0(\5f 2 ) . (53) 
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This gives two points in each of the Figures El and 0J at which the curve of vanishing twist 
emanating from the origin crosses the boundary of the image of the energy-momentum 
map. The value of the rotation number at these points is 

V5W ± = iyX± 1(5 - 2b) + 0(HV2) . (54) 

We have now treated both limiting cases, that near the focus-focus point, and that 
near the elliptic relative equilibria. The curve of vanishing twist for all values in between 
can in principle be computed from the derivative of (|36j) . and will connect the results from 
the two limiting cases for small values of 5. 



A Expansion of W 

We need to expand the Legendre standard integrals K, E, and II in the limit k — ► 1. 
For the expansion of the integral of the third kind it is important to take the inequality 
k 2 < n < 1 (|4ip into account. In this limit the following formulas can e.g. be found in pQ: 

U(n,k) = K(k) + ^irK(l- A (9,k)), (55) 



f n \ ^ f (^neg ^max)^max N \ ^ (56) 

V(l - n)(n - k 2 )) V z neg z m - m J 

Here Ao is Heumann's Lambda function, which can be expressed in terms of incomplete 
elliptic integrals F(9, k) and E(9, k) of the first and the second kind, respectively, as 

A o (0, k) = ~ (K(k)E(e, k') - (K(k) - E(k))F(9, k')) , (57) 



i 1 W . / ^min(^max -Zneg) /m\ 

arcsm \ j = arcsm \ ^ . (58) 

J- y 2i max ^.z mm z ne g) 

Note that the complementary modulus k! = y/l — k 2 and the parameter satisfy 

k' 2 = 1 - k 2 = 0(e), l-n = 0(e), 

see (|4Up. Accordingly k' 2 and 1 — n are of order e, while 9 is not small but of order 1. 
The prefactor 1Z in (|55j) cancels with the prefactor of II in (|H6|) . up to a factor of 1/2. 
Therefore we find the (still exact) formula 

W = c 1 K + c 2 E+^(l-A (9,k)) (59) 

where with ^/z max — z ncg = a/V2D + 0(e) it follows 

nci = u> + 2Bz +2C 3 +j/(2 z^) = u + G(£)) (6Q) 

\/2Z)^/z max ^neg ^ 

2B -J z max — z ne g Ba . . /oi\ 

7tc 2 = — 7 = = — + 0(e) . (61) 

V2D D 
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The elliptic integrals in the limit k' — > have a logarithmic divergence of leading order 
A = logl. (62) 

The convergent expansions in this limit are 

K{k) = A+^lk' 2 + 0{Ak' 4 ), (63) 



E(*0 = i + 5 (a - k' 2 + 0(Ak% 



(64) 



The incomplete elliptic integrals of modulus k! have regular expansions since k' — * 0, so 
that 

W W) - E { e, f ) _ / pgigf _ (| - *») * + 0<*"), (65) 

J V 1 - k /2 sm 2 w V 2 4 / 
o v 



o 



For the Heumann Lambda function this gives 



A o (0, fc) = - I 6 + - I A - - 1 sm{28)k' z + 0(Ak'*)j . (67) 

The leading order terms in the expansion of (|36|) come from the diverging K. Since E — ► 1 
for k — * 1 there is only a constant contribution from c%. From Ao the leading term is 
merely 9, so that all together 

nW = ^A + ^-9 + ^ + 0(e) 
a 2 V 

It remains to understand the parameter dependence of 6. Expanding 8 gives 



20 « 2 arcsin W - — - + 0(e). (69) 

V 2 p 

This can be simplified using the relation 



2 arcsin (3 = arctan 7 => 7 = 2— ^L^ ' 

Inserting (3 2 = {p — l)/(2p) and using p 2 = Z 2 + j 2 gives 7 = j/l, so that is proved. 
Note that when j = we have p = I and the root z m \ n collides with the pole at z = 
in the third kind integral. This is the place where the dependence on the parameters of 
n(rj, k) is continuous, but not smooth. 



16 



References 



[1] Milton Abramowitz and Irene A. Stegun, editors. Handbook of mathematical functions 
with formulas, graphs, and mathematical tables. Dover Publications Inc., New York, 
1992. Reprint of the 1972 edition. 

[2] V. I. Arnold, V. V. Kozlov, and A. I. Neishtadt. Mathematical aspects of classical and 
celestial mechanics. Springer- Verlag, Berlin, 1997. Translated from the 1985 Russian 
original by A. Iacob, Reprint of the original English edition from the series Ency- 
clopaedia of Mathematical Sciences [Dynamical systems. Ill, Encyclopaedia Math. 
Sci., 3, Springer, Berlin, 1993; MR 95d:58043a]. 

[3] Richard H. Cushman and Larry M. Bates. Global aspects of classical integrable sys- 
tems. Birkhauser Verlag, Basel, 1997. 

[4] D. del Castillo-Negrette, J.M. Greene, and P.J. Morrison. Area preserving nontwist 
maps: Periodic orbits and transition to chaos. Physica D, 91(l):l-23, 1996. 

[5] A. Delshams and R. de la Llave. KAM theory and a partial justification of Greene's 
criterion for nontwist maps. SIAM J. Math. Anal., 31(6):1235-1269 (electronic), 2000. 

[6] J. J. Duistermaat. Bifurcation of periodic solutions near equilibrium points of Hamil- 
tonian systems. In Bifurcation theory and applications (Montecatini, 1983), volume 
1057 of Lecture Notes in Math., pages 57-105. Springer, Berlin, 1984. 

[7] J. J. Duistermaat. The monodromy in the Hamiltonian Hopf bifurcation. Z. Angew. 
Math. Phys., 49(1):156-161, 1998. 

[8] H. R. Dullin and A. V. Ivanov. (Vanishing) twist in the saddle-centre and period- 
doubling bifurcation, (preprint http://arxiv.org/abs/nlin.CD/0305033), 2003. 

[9] H. R. Dullin and J. D. Meiss. Twist singularities for symplectic maps. Chaos, 13(1) :1- 
16, 2003. 

[10] H. R. Dullin, J. D. Meiss, and D. G. Sterling. Generic twistless bifurcations. Nonlin- 
earity, 13:203-224, 2000. 

[11] J. E. Howard and J. Humpherys. Nonmonotonic twist maps. Physica D, 80(3):256- 
276, 1995. 

[12] J.E. Howard and S.M. Hohs. Stochasticity and reconnection in hamiltonian systems. 
Physical Review A, 29:418, 1984. 

[13] R. Moeckel. Generic bifurcations of the twist coefficient. Ergodic Theory Dynam. 
Systems, 10(1):185-195, 1990. 

[14] J. A. Montaldi, R. M. Roberts, and I. N. Stewart. Periodic solutions near equi- 
libria of symmetric Hamiltonian systems. Philos. Trans. Roy. Soc. London Ser. A, 
325(1584):237-293, 1988. 

[15] J. Moser. Periodic orbits near an equilibrium and a theorem by Alan Weinstein. 
Comm. Pure Appl. Math., 29(6):724-747, 1976. 



17 



[16] Helmut Riissmann. Nondegeneracy in the perturbation theory of integrable dynamical 
systems. In Number theory and dynamical systems (York, 1987), volume 134 of 
London Math. Soc. Lecture Note Ser., pages 5-18. Cambridge Univ. Press, Cambridge, 
1989. 

[17] C. Simo. Invariant curves of analytic perturbed nontwist area preserving maps. Reg- 
ular & Chaotic Dynamics, 3:180-195, 1998. 

[18] A. G. Sokol'skii. On stability of an autonomous Hamiltonian system with two degrees 
of freedom under first-order resonance. Prikl. Mat. Meh., 41(l):24-33, 1977. 

[19] Jan-Cees van der Meer. The Hamiltonian Hopf bifurcation, volume 1160 of Lecture 
Notes in Mathematics. Springer- Verlag, Berlin, 1985. 

[20] Alan Weinstein. Normal modes for nonlinear Hamiltonian systems. Invent. Math., 
20:47-57, 1973. 



18 



